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ABSTRACT 

Hubble Space Telescope studies of the resolved stellar population of elliptical 
galaxies have shown that the galaxies form by steady accretion of gas which is 
all the while forming stars and evolving chemically to a metallicity distribution 
that is as high as solar composition in the most massive objects that have been 
analyzed, and much lower for low mass ellipticals. In this paper we study for 
the first time the stellar content of an early type spiral galaxy, the massive disk 
galaxy, the Sombrero, NGC 4594. We consider whether the metallicity distribu- 
tion function (MDF) in the observed field matches that of elliptical galaxies of 
some luminosity, and what these data imply for the accretion and enrichment 
model that can be fitted to the MDF. The MDF of NGC 4594 is similar to that 
of the elliptical galaxy of similar luminosity, NGC 5128. The field we are probing 
is a combination of the galaxy's bulge and halo. 

Subject headings: galaxies: abundances - galaxies: photometry - galaxies: for- 
mation 



Introduction 



Deep imaging of nearby elliptical galaxies has in r ecent years opened up their stellar 
populations for analysis of the distribution of metallicity. lHarris et al. I (120071 ) has shown that 
the metallicity distribution is measurable in a range of galaxies from the dwarf spheroidals 
to giant ellipticals such as NGC 5128, and that the the metallicity distribution functions 
(MDF) are sequenced in luminosity, such that Draco is the most metal poor in mean Z, L* 
galaxies are intermediate, and NGC 5128 the most metal rich. In massive galaxies the MDF 
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in metallicity peaks below solar Z and has an extended tail to very lo w metallicities . Lowe r 
mass galaxies show a Gaussian MDF at lower metallicities, see e.g. iLianou et al. I ( 120101 ) . 
The MDF constrains models of elliptical galaxy formation, including monolithic c ollapse 
model s, hierarchical merging models, and merger scenarios over cosmic time, see e.g. Ilkuta 



(120071 ). 



Nearby edge-on galaxies present a unique opportunity to learn about the formation 
history of the bulge, thick disk, and halo of galaxies. Seth, Dalcanton, & de Jong (2005) and 
Mould (2005) have examined HST imaging of two samples of late type edge-on galaxies and 
found thick disks without metallicity gradients and without asymptotic giant branch stars 
above the red giant branch tip. These observations rule out models of slow collapse and 
of young disk heating as thick disk formation mechanisms. Models that would be favored 
include those of Brook et al (2004) whose chemo-dynamical cosmological simulations identify 
thick disk formation with a z > 2 epoch of chaotic merging of multiple peer disks of mixed 
gas and stars. 

Bulge dominated disk galaxies can be studied in the edge-on case in the same way. A 
superb example is the Sa galaxy NGC 4594. We present such a study in this paper an d 
compare it with the Milky Way analog, NGC 891, faeikuba et al. Ibood ; llbata et al. Il2009h . 
as well as with its elliptical siblings. There is no evidence of a thick disk in the perpendicular 
profiles of iBurkhead 1 Jl986h . whose scale height is approximately 400 pc. So we concentrate 
on the bulge and halo. The Sombrero galaxy is the first early type spiral to be s o studied, 
although, in a not-so-edge-on cas e, we note the w o rk of iRichardson et al. I (120091 ) on M31, 
which is limited to the outer halo. iMouhcine et al. I ( 120051 ) has studied the disk galaxies NGC 
55, 247, 253, 300, 3031, 4244, 4258, and 4945. They find that the MDFs are dominated by 
metal- rich populations, with a tail extending toward the metal poor end and that, to first 
order, the MDFs are similar to what is predicted by simple, single- component model of 
chemical evolution with the effective yields increasing with galaxy luminosity. 



2. Photometry 



The HST archive contains six undithered F606W WFPC2 exposures of a field in the 
Sombrero galaxy and six exposures in F814W, amounting to 11,700 seconds in each filter. 
These were obtained as part of the medium deep survey of the original HST Key Projects, 
GO program 5369. The field is 6.8 effe ctive radii from the center of the galaxy and is located 
as depicted in Figure 1. According to lYoung I (119761 ) ~90% of the luminosity of an r to the 
quarter law spherical bulge is contained within this radius, and so we can expect a significant 
fraction of the stars in this field to be halo stars. 
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The frames were combined with IRAF and processed with DAOPHOT (IStetson 1119871 ). 
A sample of the field is shown in Figure 2, allowing an impression of the level of crowding. 
Point spread functions were borrowed from the extragalactic distance scale Key Project 
( IKennicutt. Mould, fc Freedman Ill995[ ) and fitted to the images with ALLSTAR. The data 
underwent a second cycle with stars found in the residual image. The ALLSTAR files were 
then combined with DAOMASTER. The number of stars detected on the four chips was 
3647, 10,030, 10,158 and 11,852, where chip 1 is the planetary camera. Nine bright stars 
w ere used t o dete rmine aperture corrections to an accuracy of 0.01 mag. The calibration 
of iDolphin I (l2000[ ) was employed to transform the color magnitude diagram (CMD) to the 
VI system. There are spatial density gradients across the WFPC2 field in the point sources 
detected of 10-20%. Both the PC and WF chips are used. The stellar density increases in 
the field in the direction towards the center of the galaxy. 

The CMD is shown in Figure 3. Stars have been excluded if their color uncertainties 
exceed 0.5 mag. The CMD shows the tip of the red giant branch (TRGB) of NGC 4594. 
The TRGB is at I = 26.2 ± 0.2 mag, corresponding to (m — M)o = 30.15 ± 0.2 mag (see 



Mould fc Sakai I ( 120081 )). The population is an old one. 



NGC 4594 is at a distance o f 29.95 ± 0.18 mag according to surface brightness fluctuation 

measurements by iTonrv et al. ( 2001 ). This is the distance we adopt. The reddening is 

E(V-I) = 0.071 from ISchlegel et al. I ( 119981 ). From the distribution of giant branch colors a 



magnitude or so fai nter than Mj = -4.0, it is thus possible to measure the MDF, following 
Harris et al. I ( 120071 ). who provide template giant branches for this purpose. The MDF is 
obtained by interpolation directly in color-magnitude space and is shown in Figure 4. The 
isochrone age was 12 Gyrs and the isochrone spacing was 0.15 dex in [Fe/H]. 

Eight hundred artificial stars with V-I = 1.75 mag and 26 < I < 28 were added to 
the F606W and F814W frames, and the measurement process was repeated. The average 
color error was 0.42 mag with little dependence on magnitude. The average I magnitude 
error was 0.44 mag. The color bias due to blending was measured at -0.02 ± 0.03, which is 
negligible. The process was repeated for V-I = 3, and the color error rose to 0.54 mag. Since 
the sensitivity of color to metallicity increases in the red, this increase is not an issue for the 
MDF. Incompleteness affects the MDF at the metal rich end. However, this effect is modest: 
a factor of two incompleteness in the solar metallicity bin in Figure 4. But a small high 
metallicity tail with [Fe/H] > may have been erased by incompleteness. Figure 5 results 
from photometry of a further 2400 artificial stars and shows completeness as a function of I 
magnitude. 
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Chemical evolution 



Since the time of iPagel Patchett I (119751 ) the MDF has been used to constrain the 
chemical evolution of galaxies, starting with the Milky Way. With deviations from the 
closed box chemical enrichment assumption it is possible to fit a metal rich distribution like 
that of the solar neighborhood with dA/dS — > a, where A is the mass of accreted gas, 
S is the mass of stars made, and a is the fra ction of low mass long-lived stars formed in a 
generation of star formation. lHartwick I (119761 ) showed that a metal poor MDF, such as that 
of Galactic halo globular clusters, could be formed with gas loss from the system, all with 
standard values of the yield, p, the fraction of heavy elements, Z created in a generation of 
star formation and return to the interstellar medium. His equation (7) can be differentiated 
to yield: 



dS_ 
dZ 



oc e 



-(a+c)Z/p 



(1) 



where c is the ratio of gas lost to stars formed, when star formation is complete. The 
resulting MDFs for different values of c are shown in Figure 6. 

We also consider a somewhat more general case in which the M DF can be used to 
constr ain A(t). If \x is the gas fraction at time t, we have equation (34) of IPagel fc Patchett 
(119751 ) (mass conservation). 



fi{t) 



aS(t) + A(t) 



(2) 



If Z' is the Z of the accreted gas, heavy element accounting (equation (6) of lPagel fc Patchett 



( 119751 )) gives us: 



fxdZ = apdS-(Z - Z')dA (3) 

Then the MDF can be written as: 

dS a , , 

— = (4) 

dZ ap-{Z - Z')dA/dS y 1 

So, given the star formation history, S(t), the accretion history, A(t), and the Schmidt- 
Kennicutt law, dS/dt oc fi n , 1 < n < 2, one can predict the MDF and compare with 
observations. An alternative way to proceed is to use Z(t) as a proxy for time and parame- 
terize the gas depletion history. 
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fi = e~ zl{px) (5) 

where x is a free parameter controlling the rate of star formation decline. We can then use 
the MDF directly to calculate the accretion rate. 

{Z - Z ' ) ^ = ap %- fi (6) 

Figure 7 shows the solution for the Sombrero's MDF. The accreted gas metallicity was 
take n to be Z' = 0.0005 and x was put equal to 0.5. Values of the yield and a were taken 



from iPagel fc Patchett I ( 119751 ). The interesting results are that there is an early burst of 
accretion and that dZ/dt is fairly steady, consistent with the assumption that Z is a good 
proxy for time. If dZ/dt is ~0. 02/13 Gyr, the duration of the initial accretion burst is ~0.17 
Gyr. Other choices of parameters yield similar results, although dZ / dt becomes monotonic 
decreasing, if the gas depletion history is rapid. 



Discussion 



Our measurement of the Sombrero's MDF is, however, severely affected by photometric 
errors. As discussed in §3, these range from 0.42 to 0.54 mag in color, moving from blue to 
red on the giant branch. This corresponds to -1, +0.3 dex at [Fe/H] = -1 and ±0.2 dex at 
[Fe/H] = -0.4. To examine this, we ran a simulation with 800 artificial stars with I = 26.4 
and V-I = 1.75. That is a delta function in metallicity with [Fe/H] = -0.875. The simulated 
MDF (generated in exactly the same way as for the real data) showed a strong low metallicity 
wing but fewer solar metallicity stars than Figure 4. We therefore ran the metallicity delta 
function through our dA/dZ calculation. The result is very similar to Figure 7 with a similar 
early accretion peak. 

It is interesting to note that the pea k of the MDF i n Fig ure 4 coincides with the metal- 
rich peak of the globular cluster MDF of lHarris et al. I ( 120101 ) . 



Harris et al. I (120071 ) have shown that the MDFs of the elliptical galaxies they have 
studied (see their figure 15) form a sequence in increasing luminosity and increasing mean 
metallicity. With a luminosity M B = -21.2 and modal metallicity [Fe /H] = -0.5, NGC 4594 
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Harris et al. 


2009) 



to it s equivalent elliptical galaxy and much more metal rich than the Milky Way analog NGC 
891 (IReikuba et al. 1120091 ) ([Fe/H] « -1.0). This suggests that similar processes are occurring 
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to determine the chemical composition of early type spiral galaxy halos and elliptical galaxies. 

Our models of chemical evolution directed at the MDFs of galaxies are consistent with 
the idea that accretion history and gas loss may be key determinants. Low mass galaxies 
are unable to retain their gas and display low mean metallicity. High mass galaxies accrete 
substantial gas early in their evolu t ion a nd process that gas to a high mean metallicity. 



The mass sequence of lHarris et al. I (120071 ) 's MDFs from the Draco dwarf, through the Leo 



ellipticals to NGC 5128 may thus be explained through the agency of gas loss and accretion. 

We have glibly spoken of the field we have observed as characteristic of NGC 4594. A 
full study of the stellar population of the Sombrero would search for chemical gradients in the 
galaxy, characteristic of chemical enrichment during collapse and seek distinctions between 
the bulge, disk, and globular cluster MDFs. We leave this for further work. 



5. Conclusions 

The MDF of a field 6.8 effective radii from the center of the Sa galaxy NGC 4594 has 
a modal metallicity of [Fe/H] = -0.5, similar to that of NGC 5128, an elliptical galaxy of 
similar luminosity. This is the first early type disk galaxy to have its MDF measured. 

Such a metal rich MDF can be achieved in accreting chemical evolution models with 
an initial burst of gas accretion lasting of order 0.17 Gyrs. This is the time elapsed between 
redshifts as high as 14 and 10. Alternatively, a loss of two thirds of the gas from this field 
during star formation, say, on to the disk could have yielded the observed MDF. 

Higher signal to noise CMDs of the Sombrero over multiple fields will help to unravel 
the star formation and chemical enrichment history of galaxy bulges. 

We thank the referee for comments which have improved the paper. This work is based 
on data obtained with the NASA/ESA Hubble Space Telescope. In addition to DAOPHOT, 
this research has made use of IRAF. IRAF is distributed by NOAO which is operated by 
AURA for NSF. 
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Fig. 1 — The field of NGC 4594 with the WFPC2 footprint shown. 
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Fig. 2. — An enlargement of a 300 x 300 pixel area in the N4594 field from the center of 
chip 2. 
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Fig. 3. — Color magni tude diagram of NGC 4594. The solid curves are the giant branches of 
M15 and 47 Tuc from lDa Costa fc Armandroff I (jl990j) moved to the distance and reddening 
of the galaxy. The dashed lines are the 50% completeness lines V = 29.3 mag and I = 27.7 
mag. M15 has [Fe/H] = -2.17 and 47 Tuc has [Fe/H] = -0.71. 
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Fig. 4. — Metallicity distribution of the detected red giants in NGC 4594. 
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Fig. 5. — Photometric completeness as a function of I magnitude. The 50% completeness 
level shown in Figure 3 is at I = 27.7 ± 0.3 mag. 
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Fig. 6. — Metallicity distribution in the simple model of chemical evolution with gas loss. 
The ratio of gas lost to stars formed takes the values for distributions from left to right of 
32, 16, 8, 4 ,2, 0. 
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Fig. 7. — The accretion rate as a function of metallicity (top left) in NGC 4594. The 
independent variable, [Fe/H] = \ogZ/Z &1 can be considered to be a proxy for time. The 
assumed gas depletion rate is depicted in the top right. The chemical enrichment rate 
(arbitrary units) is at the bottom left, and the star formation rate (arbitrary units) at the 
bottom right. 



